Fidelity in topological superconductors with end Majorana fermions 
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Fidelity and fidelity susceptibility are introduced to investigate the topological superconductors 
with end Majorana fermions. A general formalism is established to calculate the fidelity and fidelity 
susceptibility by solving Bogoliubov-de Gennes equations. It is shown that the fidelity susceptibility 
manifest itself as a peak at the topological quantum phase transition point for homogeneous Kitaev 
wire, thus serves as a valid indicator for the topological quantum phase transition which signals the 
appearance of Majorana fermions. The effect of disorders is investigated within this formalism. We 
consider three disordered systems and observe fidelity susceptibility peak in all of them. By analyzing 
the susceptibility peak, we notice that the local potential disorders and the hopping disorders can 
shift the phase transition point, while off-diagonal disorders have no obvious influence. Our results 
confirm that the existence of topological quantum phase transition is robust to these disorders, while 
the behavior of the phase transition might be influenced by disorders. 

PACS numbers: 74.25.Dw, 03.67.-a, 74.62. En 



I. INTRODUCTION 

Majorana fermions (MFs), fermions vifhich are their 
own anti-particles, are under intense investigation re- 
cently. MFs have been predicted as zero-energy bound 
states in spinless p-wave superconductorsii^, known as 
topological superconductors. In these topological su- 
perconductors, two MFs form one topological qubit 
which resists to local electro-magnetic disturbance, thus 
might be useful for long-time storage of quantum 
informationii^i^. Moreover, braiding of these MFs rotates 
the topological qubits they form^'^'^, thus may constitute 
as one of the cornerstones of quantum computatiori^i^iir— . 

MFs were originally proposed in pure theoretical 
toy modelsii^. However, recent progress in spin-orbit 
coupling research makes it possible to be realized in 
hybrid systemaiS, such as superconductor-topological 
insulator interface^i, or semiconductor-superconductor 
heterostructur6i2r2&. Among these systems, one- 
dimensional spin-orbit coupling nanowire in proximity 
to conventional superconductor is of particular interest, 
due to its relatively simple structurei^ii^. It was pre- 
dicted that there is a topological quantum phase transi- 
tion (TQPT) in the system when proper Zeeman field is 
applied, and zero energy end MFs appear in the topo- 
logical non-trial phase. Recently, a device with required 
condition has been fabricatedi^, and the differential con- 
ductance was measured under various strength of Zeeman 
splitting. A robust zero-bias peak has been reported, in 
agreement with the resonant Andreev effect which was 
proposed to signify the MFs^^. Moreover, the Zeeman 
energy in the device was manipulated through rotating 
the applied magnetic field, and the zero-bias peak was 
discovered under the exact condition predicted by the- 
ory. Later on, subsequent experiments have been con- 



ducted on similar devices, and their results validated the 
previous repor t^^i^° . 

The realization of MFs seems very optimistic with 
these experimental progresses. However, there are some 
issues need to be fixed before the MFs can be completely 
pinned down. In the previous experimental reports^^, the 
extracted superconducting energy gap does not close at 
the TQPT point where the zero-bias peak appears, more- 
over, the height of the zero-bias peak is much lower than 
the quantized conductance h/2e predicted by the reso- 
nant Andreev effecti^. Because of these discrepancies, 
there are arguments on the proper interpretation to the 
recent experimental results. Some provide theoretical ex- 
planations for these two issues, and try to certify the ex- 
istence of MFs in the system^i^— . Others point out that 
the appearance of zero bias peak is due to non-topological 
mechanisms, such as disorders^ir— , edge confinement^, 
or Kondo effect^. While the second issue can be at- 
tributed to finite temperature effect, the first issue is 
rather difficult to understand. Traditionally the energy 
gap which protects the topology of system should close 
at the TQPT point, where the spectrum is reconstructed 
and the zero energy MFs may appear. This discrepancy 
indicates that the TQPT in the topological supercon- 
ductors is an important problem, and might be crucial 
to understand the zero-bias peak. Therefore, more inves- 
tigations on the TQPT are necessary to understand the 
appearance of zero bias peak and its robustness. 

One established method to investigate quantum phase 
transition is the fidelity approach^^. Fidelity is a measure 
of the difference between two quantum states. Intuitively, 
it is a natural marker for quantum phase transition, since 
it should presents a drop at the phase transition point 
where a dramatic reconstruction of the quantum states 
happens. In application, the second derivative of the fi- 
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delity, i.e. fidelity susceptibility22-, would be a better 
choice for application. It should manifest itself as a sharp 
peak at the quantum phase transition point, where the 
peak height scales with the system size. Fidelity and fi- 
delity susceptibility can describe quantum phase transi- 
tion without assuming a priori knowledge of the system, 
thus are particularly suitable for dealing with TQPT, 
where the traditional symmetry-breaking formalism does 
not work^. 

In this work, we apply the fidelity and fidelity sus- 
ceptibility to study the zero temperature TQPT in one 
dimensional Kitaev wire, which supports end Majorana 
fermions at the topological non-trivial phase. We estab- 
lish a scheme to calculate the fidelity and fidelity sus- 
ceptibility by solving the Bogoliubov-de Gennes (BdG) 
equations. We confirm the validity of this scheme by com- 
paring the results from BdG equations in open boundary 
system and the analytic results in periodic boundary sys- 
tem. In both systems, the fidelity susceptibility manifest 
itself as a peak at the phase transition point predicted 
from the topological arguementi, thus serves as a valid 
marker for the TQPT. Then we introduce this scheme 
to investigate the disordered systems, where no informa- 
tion for the quantum phase transition has been drawn 
analytically. We consider several simple disorders, such 
as the local potential disorder, the hopping disorder, and 
the off-diagonal disorder. We find that the fidelity sus- 
ceptibility always manifest itself as a peak for all dis- 
ordered systems under consideration, implying that the 
existence of the TQPT is robust to these disorders. This 
result agrees with previous topological arguments, and 
supports the conclusion that the MPs are robust to dis- 
orders. After analyzing the position and the height of 
the peak in fidelity susceptibility, we are able to provide 
detailed information for the TQPT in disordered system 
besides its existence. It is found that the local poten- 
tial disorders and the hopping disorders shift the phase 
transition point obviously, while the off-diagonal disor- 
ders have no obvious infiuence. We also perform scal- 
ing analysis on the TQPT, and find that the quadratic 
scaling behavior for homogeneous system is preserved in 
local potential disordered system and off-diagonal disor- 
dered system, while the scaling behavior of the hopping 
disordered system significantly deviates from the simple 
quadratic form. From these results, we are able to tell 
that the hopping disorder has a strong influence on the 
TQPT, while the off-diagonal disorder has little influence 
on the TQPT. Our work illustrate that the fldelity ap- 
proach is very useful to find the TQPT in topological 
superconductors, thus can provide an evidence for the 
existence of MFs, since MPs are always accompanied by 
a TQPT which is absent for accidental zero energy exci- 
tations. For all theoretical systems where zero energy ex- 
citations appear and argued to be MFs, fidelity approach 
can be applied and provides a useful check by detecting 
the TQPT in the system. From this point of view, our 
study can help to resolve the debate over the mechanism 
behind the zero bias peaks discovered recentlyii. 



This paper is organized as follows: The model and the 
analytic solution for the periodic boundary condition are 
discussed in Sec. II. The scheme for calculating fidelity 
and fidelity susceptibility by solving BdG equations, and 
its application to open boundary condition are presented 
in Sec. III. The disordered systems are solved in Sec. IV. 
Finally, we give a summary in Sec. V. 



II. MODEL AND ANALYTIC SOLUTION 

The topological superconductor fabricated recently is 
a hybrid-systemii, where a spin-orbit coupling nanowire 
under a Zeeman field is in proximity to a conventional 
s-wave superconductor. A detailed model for this de- 
vice should involve many factors such as Zeeman energy, 
spin-orbit coupling, charging energy etc. However, its 
key physics can be captured by the one dimensional Ki- 
taev wire^, which is a one dimensional lattice model for 
a spinless superconductor, 

N N 
N 

+ 5](Ac,c,+i+A*c]^,c]), (1) 

where is the electron creation operator, t is the hop- 
ping integral, /j, is the chemical potential, A is the su- 
perconducting gap, and N is the lattice size. It has been 
established by Kitaev that a topological non-trivial state 
with end MFs appears in this model^, and the TQPT 
appears at ^ = 2i for |A| > 0, where the energy gap of 
the system closes. While most research concentrated on 
the topological non-trivial state, only few of them have 
targeted on the TQPT itself. In this paper, we apply 
the fidelity approach to investigate the TQPT in this 
model. Let us first look at the the simplest case when 
the wire formes a ring, imposing a periodic boundary 
condition to the model. The zero energy MFs vanish in 
this ring geometry, however, the two topological states 
and the TQPT between them persist. With the peri- 
odic boundary condition, the model can be simplified via 
the Fourier transformation, where momentum is a good 
quantum number. A pair of electrons with momentum 
k and —k will form the Cooper pair, and a conventional 
BCS ground state wave-function can be obtained. 
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where quasi-particle coefficients are, 
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FIG. 1: (Color online). Fidelity (a) and Fidelity susceptibility 
(b) as a function of chemical potential /i, and (c) the height 
of peak in fidelity susceptibility as a function of lattice size, 
where A = t. 



From this energy spectrum, it can be seen that the energy 
gap closes at = 2t, if |A| > 0. This special point where 
the energy gap vanishes is argued to be the TQPT point 
of the system^, since topology of the system is protected 
by the energy gap. 

Choosing the chemical potential /i as the driving pa- 
rameter, the fidelity between two nearby ground states is 
obtained straightforwardly with this wave-function, 

/ ={^G{f^)\^G{p^ + Spi)) (5) 

= Yl_ [ufe(M)ufe(M + (^m) + Vfc(^)vfc(/i + Sn)] , 
fc 

and the fidelity susceptibility is the derivative of the fi- 
delity, with an analytical form. 



XF 



2 lim 

sin^ k 



l-|(vI/G(/i)|*G(M + 5M))l' 



,d*G 



' d^l 



(6) 



The fidelity susceptibility is simply determined by the 
energy spectrum, however, this simple expression is ap- 
plicable only under the periodic boundary condition. 

The results of fidelity and fidelity susceptibility with 
a system size N = 1000 are shown in Fig. la. and Fig. 



FIG. 2; (Color online). Same as in Fig. 1 except for open 
boundary condition. 



lb. as a function of chemical potential /z. It can be 
seen that the fidelity susceptibility manifests itself as a 
sharp peak at /i = 2i, and clearly signatures the TQPT 
point speculated in previous researchesi. Our results con- 
firm the previous argument based on the energy spectrum 
analysis^, and validate the fidelity approach for marking 
the TQPT in this system. 

We also perform a scaling analysis for the fidelity sus- 
ceptibility at the TQPT point. The height of the peak 
is shown as a function of the lattice size N in Fig. Ic. 
It is obvious that the height of the peak quadratically 
scales with the system size, and would diverge when ap- 
proaching the thermodynamic limit. We notice that our 
result for Kitacv model should be the same to the previ- 
ous results obtained in the transverse-field Ising model^^, 
since these two models are equivalent to each other in the 
absence of disorders. 



III. OPEN BOUNDARY PROBLEM 

The most fascinating aspect of topological super- 
conductor is the topologically protected end Majorana 
fermions, which unfortunately is missing in the ring-like 
structure as we discussed in the previous section. To get 
the zero energy MFs, we must consider a wire with ends 
which enforce an open boundary condition. In this case, 
no analytical solution is available, due to the absence of 
the translational symmetry. A conventional method to 
tackle this system is to solve the Bogoliubov-de Gennes 
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(BdG) equations, diagonalize the Hamiltonian in real 
space and obtain the energy spectrum as well as the 
quasi-particle wave functions. In order to calculate the fi- 
delity and fidelity susceptibility, we have to construct the 
wave function of the system. However, a wave function in 
the electron representation like the one given in Eq. (2) is 
absent, since the pairing does not happen in momentum 
space. On the other hand, the ground state wave func- 
tion can be constructed within the quasi-particle repre- 
sentation simply by filling all the negative energy quasi- 
particle states, similar to the picture of fermi sea. 

In the following, we utilize this formalism and find that 
the fidelity and fidelity susceptibility naturally arises. 
The system we study here is a finite size Kitaev wire with 
N sites. The energy spectrum En and the quasi-particle 
wave- function ipnij) = (unj, v„j) can be obtained by nu- 
merically solving the BdG equations, where n runs from I 
to 2N denoting the quantum number of the energy level, 
and j runs from I to A'^ denoting the site number. As we 
have discussed, the ground state wave function should 



be the combination of the wave functions of all negative 
energy quasiparticle states, which can be expressed as a 
Slater determinate. 



l*0-l,---,Jn)) = 



V'n(jl) ••■ V'n(in) 



(7) 



where n denotes the number of quasi-particle states with 
negative energy. For the open boundary Kiatev wire, 
there are zero energy MF bound states at both ends in 
topological non-trivial phase, making the ground state 
doubly degenerate. In order to avoid complexity, we 
added quasi-particle wave function of these zero energy 
states into Eq. (7), making n = A'^ -I- 1 in our calculation. 

Inserting this wave function into the definition of fi- 
delity Eq. (5), and taking advantage of a simplification 
for calculating fidelity in multi-particle fermion system. 
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where indicates the wave-function under chemical po- 
tential fi, and Pim = J dx'il;*^i{x)iljf_i+Sfi,m{x). With this 
expression for fidelity, the fidelity susceptibility can be 
obtained through the definition 
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where the differentiation is calculated numerically. The 
result of the fidelity and fidelity susceptibility are shown 
in Fig. 2a and Fig. 2b, and the scaling behavior for the 
peak of the fidelity susceptibility is shown in Fig. 2c, 
with all physical parameters the same as in the previous 
section. We find that the fidelity susceptibility manifest 
itself as a sharp peak at the TQPT point, showing that it 
is capable of signifying the TQPT in the wire. The height 
of the peak is quadratically scaling with the system size, 
similar to the results discovered in Fig. Ic which comes 
from the analytic solutions in periodic boundary condi- 
tion. Considering that the boundary condition should be 
irrelevant to the phase transition for large system, this 
near identical result illustrates that our scheme of cal- 
culating fidelity by solving the BdG equations is a valid 
procedure, and would be helpful to investigate the TQPT 
in the superconducting systems. 



IV. DISORDER EFFECTS 

After establishing the fidelity approach to the TQPT 
by solving the BdG equations, we are now ready to deal 
with more realistic and complicated disorder systems. 
The role of disorders in experimental device is still under 
debate. In traditional believe, its effect should be very 
weak. However some recent theories argue that disorder 
itself might be the origin of the zero-bias peak, making 
the zero-bias peak irrelevant to the MFs^i"— . In this 
section, we investigate the effects of the disorders on the 
TQPT based on the scheme we built in last section. Since 
TQPT is necessary for the appearance of MFs, our study 
will also provide the information on the robustness of the 
MFs. 

We study some examples to show the power of fidelity 
approach to investigate the TQPT in disordered system, 
and illustrate the influence of several typical disorders on 
the TQPT. The most typical disorder is a local potential 
disorder, which modulates the local chemical potential 
with the Hamiltonian, 
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FIG. 3: (Color online). Fidelity (a) and Fidelity (a) suscep- 
tibility as a function of chemical potential /i with disorder 
concentration of 20% (solid line), 50% (dashed line), and 80% 
(dotted line), where lattice size = 200. (c) The fidelity sus- 
ceptibility of with lattice size — 100, 200, and 400, where 
the disorder concentration is 50%. (d) The height of the peak 
in fidelity susceptibility as a function of the lattice size A^, 
with 50% of disorders. 



where i indicates the site with the impurity, and Vi in- 
dicates the local potential. We consider a random dis- 
tribution of the same kind of disorders on the wire, and 
investigate the behavior of the TQPT for different dis- 
order concentrations. The two mathematical limits of 
one hundred percent of disorders and zero percent of dis- 
orders are trivial and well understood theoretically. In 
the first limit, an extra potential should be enforced to 
every site, which is equivalent to shifting the total chem- 
ical potential from /i to fj,— Vi, thus the phase transition 
point should be shifted from fi — 2t to jj. = 2t + Vi. In 
the second limit, there is no disorder in the system, and 
the phase transition point should be the same as in the 
previous section. However, the properties of the TQPT, 
even its existence, for a intermediate concentration has 
to be investigated within fidelity approach. The results 
of the fidelity and fidelity susceptibility with 20, 50 and 
80 percent of potential disorders are shown in Fig 3a, and 
Fig. 3b, with the disorder strength Vi — t. We can see 
that for all disorder concentrations, fidelity susceptibility 
manifest itself as a peak at different places. It is a clear 
evidence that the TQPT is robust to disorder, while the 
phase transition point shifts when changing the disorder 
concentration. In order to further confirm the robustness 
of the TQPT, we present the fidelity susceptibility for dif- 
ferent lattice size of 100, 200, and 400 in Fig. 3c, with 50 



FIG. 4: (Color online). Same as in Fig. 3 except for hopping 
disorders. 



percent of disorders. It can be seen that the sharp peak 
in fidelity susceptibility appears for all lattice sizes, while 
the peaks are located at almost the same positions. We 
argue that the small deviation of peak positions might 
come from different specific disorder configurations, and 
reflect the finite size effect. We present the scaling be- 
havior of the peak height in Fig. 3d, and find it roughly 
scales quadratically with the lattice size. However, we 
have to mention that this scaling behavior is very inac- 
curate, since the peak height still depends on the specific 
disorder configurations within the lattice size we have 
investigated. We have performed more calculations for 
different lattice sizes and different disorder strengths, as 
well as ensemble averages for the disorder, and find that 
the peak in fidelity susceptibility is very robust. These 
results clearly suggest that the TQPT, as well as the 
appearance of MFs, will not be suppressed by the local 
potential disorders. 

The second type of disorder we investigate is the hop- 
ping disorder, which modulates the local hopping integral 
between sites with the Hamiltonian, 

H2 = V2clc,+i + h.c. (11) 

where the disorder is located at the bond between site 
i and i + 1, with V2 the disorder strength. Similar to 
the local potential disorders, the limiting case of the dis- 
order concentration is trivial. For one hundred percent 
of disorders, the phase transition point is simply shifted 
to ^ = 2(t — V2). We present the results of fidelity and 
fidelity susceptibility with 20, 50, and 80 percent of hop- 
ping disorders in Fig. 4a and Fig. 4b , where the disorder 
strength is V2 = 0.5t. We can see that the influence of 
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FIG. 5: (Color online). Same as Fig. 3 except for ofF-diagonal 
disorders. 



strength V-^ = 0.5t. As we expected, we find that the fi- 
delity susceptibility manifest a sharp peak almost at the 
same position as in the homogeneous system, indicating 
a minimal influence to the TQPT. We also present the 
fidelity susceptibility for different system sizes in Fig. 5c 
with 50 percent of disorder, and find that the peak al- 
ways appears. The scaling behavior of the peak height is 
shown in Fig. 5d, and it is found that the height of the 
peak is scaling quadratically with the system size, which 
is very close to the result for homogeneous system in the 
last section. These results illustrate a minimal influence 
of these off-diagonal disorders, and confirmed that they 
are irrelevant to the TQPT and the MFs. 

We have studied three types of the disorders in this 
section and found that all of them preserve the TQPT. 
While the on-site potential disorder and hopping disorder 
shift the phase transition point and modulate the scal- 
ing behavior of peaks in fidelity susceptibility, the off- 
diagonal disorder has a minimal infiuence to the TQPT. 
Our results show that the TQPT, together with the topo- 
logical non-trivial state and the end MFs, is robust to 
these disorders. 



hopping disorder is very similar to the local potential 
disorder, that is. The TQPT is robust while the phase 
transition point is shifted. The fidelity susceptibility for 
different system size is shown in Fig. 4c, and it is found 
that the peak always appears indicating the existence of 
a TQPT. The scaling behavior of the peak height is pre- 
sented in Fig. 4d, and it is found that the quadratic 
scaling for homogeneous system barely exists anymore. 
These results illustrate that the hopping disorder pre- 
serves the TQPT in the system, however, it has a larger 
influence on the detailed behavior of the TQPT than the 
local potential disorders, which can be seen from the scal- 
ing behavior of height of the peak in fidelity susceptibility. 

Finally, we are going to investigate the off-diagonal dis- 
order, which has been extensively studied in high temper- 
ature superconductors^. The off-diagonal disorder can 
be understood as a static fluctuation of the amplitude of 
superconducting gap. In a mean-field level, this disorder 
can be described with the Hamiltonian, 

Hz - -yAA^i + h.c. (12) 

where the disorder is located at the bond between sites 
i and i -I- I, with V3 the modulation to the local super- 
conducting gap amplitude. For this off-diagonal disorder, 
the two limit of the disorder concentration is very similar, 
since the TQPT should be irrelevant to detailed super- 
conducting gap values as long as it is non-zeroi. In this 
case, we would expect that the TQPT is less influenced 
by these disorders. In Fig. 5a and Fig. 5b, we show the 
results of fidelity and fidelity susceptibility for 20, 50, 
and 80 percent of off-diagonal disorders, with a disorder 



V. SUMMARY AND DISCUSSION 

In summary, we have investigated the topological 
quantum phase transition in the Kitaev model within the 
fidelity approach. A formalism to calculate the fidelity 
and fidelity susceptibility by solving the BdG equations 
has been established and its validity is confirmed by com- 
paring to the analytic results in homogeneous system. 
This formalism is then applied to investigate the disor- 
dered systems, where no analytic solution is available. 
Three typical disordered systems are investigated, and 
it is found that the disorders preserve the TQPT based 
on the appearance of the peak in fidelity susceptibility. 
From the peak positions, we find that the diagonal dis- 
orders shift the phase transition point clearly, while the 
off-diagonal disorders has no obvious impact on the peak 
position. 
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